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ABSTRACT 


The utility of hyperspectral imagers for the identification, classification and status 
of a specific material based on it's spectral characteristics has been demonstrated in the 
fields of geology, forestry and meteorology. The United States military has an interest in 
the utility of hyperspectral imagers for a multitude of tactical and strategic purposes. The 
Hyperspectral MASINT Support to Military Operations Program (HYMSMO) was 
designed to explore this arena through a series of planned collection operations utilizing 
the Hyperspectral Digital Imagery Collection Experiment (HYDICE) sensor in addition to 
other hyperspectral imaging platforms. Operation DESERT RADIANCE explored many 
of the areas relating to the tactical detection and classification of military targets. 
Through the use of the Low Probability of Detection and the Principle Components 
Transformation algorithms contained in the HYDICE Starter Kit and ENVI software 
package, this thesis shows that the detection of a tactical target by use of it’s unique 
spectral signature is feasible. 
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I. INTRODUCTION 


The field of imaging spectrometry, also known as hyperspectral imaging, has 
slowly emerged over the past two decades. Hyperspectral imaging refers to remote 
sensing based collection of temporally-registered imagery over a large number of 
discrete, contiguous bands whereby a complete reflectance spectrum is obtained. With 
parallel advances in optics, microelectronics and image processing software, 
hyperspectral imaging has matured and can provide detailed information concerning the 
earth's environment and the warfighter's battlespace. 

Electromagnetic energy is the most abundant energy form available for remote 
sensing applications. Within the electromagnetic spectrum, the visual and solar reflected 
infrared regions can reveal specific reflective and emissive spectral signatures of the 
imaged object. The value of hyperspectral imaging systems is their ability to acquire 
complete reflectance spectrum for each picture element (pixel) in the image. Many 
materials have definitive absorption features that are 20-40 nanometers wide and thus are 
undetectable by the current broadband multispectral sensor systems. Hyperspectral 
sensors nominally acquire data in 10 nm bands providing sufficient spectral resolution for 
direct identification of those materials. 

Military applications for hyperspectral sensors are just beginning to be evaluated. 

In the past, manual analysis of stereo images by image interpreters has been used for 
determining the shape, size and arrangement of specific targets and their features. 
Examples include targets such as roads, airports, dams, vehicles, structures and urban 
areas. The location, identification and assessment of military significant targets within 
the battlefield can be enhanced with the use of hyperspectral sensor systems. By 
applying hyperspectral imaging techniques and analysis, other differences, not based on 
shape and size, such as color, spectral reflectance and spectral emittance (temperature) 
can be identified within the image. Specific camouflage types, including paint and fabric, 
soil composition, downed pilot taggants, stressed vegetation, thermal plumes, nuclear, 
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biological and chemical agents are a few among a list of many possible candidates for 
detection/identification by hyperspectral sensor systems. 

The data obtained for military use by hyperspectral sensors will most likely fall 
within the field of measurement and signature intelligence (MASINT). Imagery 
interpretation (II) and imagery intelligence (IMINT) products have been the traditional 
methods for evaluating and disseminating information relating to the description of 
objects, activities and terrain within the battlespace. Literal exploitation of imagery deals 
with the extraction of information from imagery by a human interpreter. MASINT, on 
the other hand, deals with the nonliteral exploitation of imagery. A hybrid definition of 
nonliteral imagery exploitation can be found in Joint Pub 1-02: 

The process of extracting non-spatial information from image data, 
automatically or semi-automatically, using non-traditional, advanced 
processing techniques, employing models, measurements, signatures 
(spectral, textural, temporal, polarization), or other features to detect, 
locate, classify, discriminate, characterize, identify (material, unit, 
function), quantify (material, time, physical), track, predict, target, or 
assess objects, emissions, activities, of events represented in the imagery. 

The Hyperspectral MASINT Support to Military Operations (HYMSMO) program 
grew out of military requirements for Operational Measurement and Signature 
Intelligence (OpMASINT) support and is sponsored by the Central MASINT Technology 
Coordination Office (CMTO). It is designed to demonstrate and advocate the utility of 
imaging hyperspectral spectroradiometers as true MASINT resources (as opposed to 
conventional IMINT assets). In addition, the HYMSMO program will examine the 
nontraditional, time-sensitive, MASINT application termed support to military operations 
(SMO). 

The Hyperspectral Digital Imagery Collection Experiment (HYDICE) is a 
Congressionally funded "dual-use" initiative that has developed an airborne, imaging 
hyperspectral spectroradiometer, offering high spectral fidelity, spatial resolution, and 
performance. HYDICE, as the sensor is called, is an integral part of the HYMSMO 
program and will be utilized for data collection experiments in the spring of 1995. Data 
2 



collection experiments using non-HYDICE sensors have been completed and will 
continue throughout the program lifetime. The purpose of these experiments is to begin 
to evaluate the use of hyperspectral sensors in the detection and classification of military 
significant targets. One such experiment. Operation DESERT RADIANCE, took place at 
the White Sands Missile Test Range in New Mexico. The purpose of this experiment 
was to demonstrate the potential for using hyperspectral remote sensing to differentiate 
among natural and man-made materials. 

The thesis is broken down into five chapters and two appendices. Chapter I, the 
Introduction is followed by Chapter II, Background, which addresses the history and 
accomplishments of hyperspectral imagery along with a brief explanation of the most 
current hyperspectral image processing and analysis algorithms dealing with material 
detection. Chapter III, Experiment Description, discusses the Operation Desert Radiance 
experiment where the data analyzed in this thesis was collected. Chapter IV, Experiment 
Observations, highlights significant observations made during the author's participation in 
the experiment. Chapter V, Data Analysis, describes the methods used in the detection of 
an actual tactical target by use of it's spectral signature and the results attained. Chapter 
VI, Conclusions is followed by the Appendix which includes all figures identified in the 
thesis. 
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II. BACKGROUND 


Hyperspectral imaging is one of the many forms of remote sensing available today 
for extracting information relating to the earth and it's surroundings. Because all 
materials reflect, absorb, or emit photons in ways characteristic of their molecular 
makeup, a high resolution trace of the intensity of the transmitted, reflected, emitted or 
luminesced radiation versus wavelength forms a graphical record unique to a given 
material. (Rinker, 1990). Figure 2.1 is a plot of the percent reflectance vs. wavelength of 
two fabrics and a green leaf. The water absorption features at 1.4 and 1.9 micrometers 
are clearly visible and the plot shows a close similarity between fabric A and green leaf 
B. Fabric C is closely matched in spectral shape out to about 1.1 micrometers but at 
approximately 1.8 micrometers where it begins to deviate significantly. These 
differences provide the means to identify materials and detect man-made objects in large 
"natural" scenes. Hyperspectral imaging refers to the imaging of an area or "scene" over a 
large number of discrete, contiguous spectral bands such that a complete reflectance 
spectrum can be obtained for the area being imaged. The field of hyperspectral imaging 
has evolved over the last 2 decades primarily due to efforts by the National Aeronautics 
and Space Administration (NASA), the Department of Defense (DOD) and recently, 
private industry research and development. 

A. HISTORY OF HYPERSPECTRAL IMAGING 

1. Concept 

Terrestrial remote sensing from aircraft and space was developed to provide 
specific information about the makeup of the earth's surface which was not readily 
obtainable from direct observation on the surface. Panchromatic sensors provide 
information in parts of the visible portion of the electromagnetic (EM) spectrum only. 
The level of detail within the image is dependent on the sensitivity of the sensor's 
detectors to differences in signal strength as it records the radiant flux reflected or emitted 
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from the object being sensed. For example, the French SPOT (System Pour l'Observation 
de le Terre) Panchromatic sensor uses 6 bit. Differential Pulse Code Modulation (DPCM) 
which provides the capability to represent the image scene in 64 different levels of gray. 
While this provides adequate detail for most human analysis techniques, it focuses only 
in the visible portion of the EM spectrum and does not exploit the differences in the 
reflectivity, as a function of wavelength, that can be represented by multispectral sensors. 

The LANDSAT Thematic Mapper (TM) is a multispectral sensor capable of 
imaging in 7 bands of the visible, near infrared (NIR), shortwave infrared (SWIR) and the 
thermal infrared portions of the EM spectrum. Figure 2.2 shows the location, in 
wavelength, of the TM bands. Note the coarse resolution and the gap from 0.9 to 1.5 
microns. Data front several bands can be combined to produce an image that contains 
information from each bandpass selected. A false-color composite image can be 
constructed by the assignment of one of the three primary colors to each of the selected 
bands. This image can, in many instances, provide additional information beyond that of 
a panchromatic image alone. 

The main limitations associated with multispectral sensor designs is their inability 
to acquire the information necessary to identify the composition of surface materials, 
because of the poor spectral resolution of the sensors. The spectral reflectance and 
emittance characteristics for surface materials, which are determined by electronic and 
vibrational energy states within the materials, are usually too highly structured to be 
observed at coarse spectral resolutions (Vane, 1985). Figure 2.2 shows characteristic 
spectra from various minerals, along with the bandpass ranges of the LANDSAT TM. 
Different causes of absorption can be identified in the mineral spectra. Absorption due to 
electronic transitions can be detected at 0.92 micrometers. Vibrational absorption 
features due to bound and unbound water become more important in the determination of 
the target material at around 1.0 micrometers. These features may however be masked 
due to the presence of strong atmospheric water bands at 1.4 and 1.9 micrometers. Many 
materials and minerals exhibit distinctive absorption features in the region from 2.1 to 2.4 




micrometers. This is due primarily to the combination bending-stretching overtones of 
the OH vibration. Laboratory work in the field of spectroscopy demonstrated the wealth 
of information contained in the reflectance spectra of various minerals and vegetation 
types thereby motivating the development of imaging spectrometry (Vane and Goetz, 
1993). A graphical depiction of the hyperspectral imaging concept is shown in Figure 
2.3. Images taken simultaneously in 100-200 spectral bands form an image cube, with 
each layer of the cube representing a different wavelength. Each pixel will have an 
continuous spectrum that can be used to uniquely identify surface materials. 

2. Evolution 

Laboratory and field experiments conducted in the 1960's and 1970's confirmed the 
possibility that solar reflectance measurements by airborne sensors could permit the 
identification of the mineral makeup of rocks, soils and vegetation. However, the spectral 
resolution for specific identification would have to be on the order of 10 to 20 nanometers 
covering the spectral region of 400 to 2500 nanometers. These results motivated the 
development of airborne and spacebome hyperspectral imagers in the late 1970's. 
Curiously, one of the first systems flew on an early Shuttle flight. 

The Shuttle Multispectral Infrared Radiometer (SMIRR) was flown on STS-2 in 
1981. SMIRR measured spectral reflectance in 10 spectral bands with a spatial resolution 
of 100 meters. Investigation of the data by A.F.H. Goetz and L.C. Rowan concluded that 
the SMIRR had detected a previously unmapped mineral deposit containing specific clay 
and carbonite minerals (Taranik et al, 1993). 

In 1980, the Jet Propulsion Laboratory (JPL) began the development of the Thermal 
Infrared Multispectral Scanner (TIMS). Imaging in the 8 to 12 micrometer range with a 
spatial resolution of 20-30 meters, depending on flight altitude, TIMS first flew over 
Cuprite, Nevada gathering thermal emission spectra of the various mineral deposits. In 
August 1994 TIMS was flown on board a NASA C-130 aircraft over the McDowell 
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Mountains and urban areas of Scottsdale, Arizona to collect data that will be used by the 
city in the development of it's Geographic Information System (GIS) (Constance, 1994). 

a. Airborne Imaging Spectrometer (AISI/II) 

JPL began development of the AIS in the early 1980's after the successful flight 
of SMIRR on the shuttle. A critical area of sensor development had been that of the 
detector which is the device that converts the energy of the incoming photons to an 
electrical signal. Early multispectral sensors, the LANDSAT Multispectral Scanner and 
the Thematic Mapper, incorporated discrete detector elements which were scanned across 
the earth's surface as the spacecraft moved along it's orbital path. The main limitation to 
these early designs was the short residence time of the detector in each instantaneous field 
of view (IFOV) thereby reducing the signal to noise ratio (SNR) of each detector element 
(Vane and Goetz, 1988). The use of multiple line arrays has been demonstrated by the 
French SPOT (System Probatorie Pour l'Observation de la Terre) satellite. An area array 
design was first utilized with AIS I. Designed as a testbed for the future development of 
space based, infrared imaging spectrometers, AIS I incorporated a 32 x 32 element, 
mercury cadmium telluride (HgCdTe) array. The AIS I sensor used a grating 
spectrometer to separate the incoming signal into 128 contiguous bands in the spectral 
region from 1.2 to 2.4 micrometers with a spectral resolution of 9.3 nanometers. AIS II 
was configured with a larger, 64 x 64 element, HgCdTe array and covered the region 
from 0.8 to 2.4 micrometers at a spectral resolution of 10.6 nanometers. Figure 2.4 
depicts the evolution of different image acquisition techniques. 

Kruse, 1987, utilized data sets acquired by both the AIS I and AIS II in an 
analysis of hydrothermally altered rocks in the northern Grapevine Mountains, of 
Nevada and California. By calculating the band position, band depth and bandwidth of 
the strongest absorption features of each spectrum and mapping these parameters into a 
hue, saturation, intensity (HSI) coded image, he was able to rapidly identify different 
mineral groups, as shown in Figure 2.5. Individual spectra from the images were 




compared with laboratory spectra obtained from on-site samples. Positive identification 
of sericite, montmorillonite, calcite and dolomite was accomplished using the AIS data. 
As shown in Figure 2.6, the distinct absorption features at varying bandpasses provide a 
means to uniquely identify a specific mineral. 

Some limitations in the use of the AIS II were found in a study conducted by 
Riggs and Running (1990). While studying Norway Spruce and White Pine plots, they 
determined that high spectral resolution imagers provided only limited information 
concerning water stress in conifers at the landscape scale. Figure 2.7 shows the relative 
reflectance plots of stressed and control White Pines. The plots clearly show that there is 
not a notable increase in the reflectance of the stressed group even though the group had 
been stressed for approximately 3 weeks. The weakness in the use of hyperspectral data 
for this purpose does not appear to have been sensor related. The problem appears to be 
that the apparent reflectances of the vegetation do not change appreciably with 
deprivation of soil water over extended periods. Only trees approaching lethally, low 
water content levels will be distinguishable. 

Forest canopy characteristics and nitrogen cycling were analyzed using AIS data 
acquired over a series of Wisconsin forest ecosystems. Analysis revealed that the canopy 
lignin or a closely associated property was influential in canopy reflectance within the 1.2 
to 1.6 micrometer spectral region. These results prove that further study of forest canopy 
chemistry with remote sensing should be conducted (Wessman, et al, 1987). 

b. Airborne Visible/Infrared Imaging Spectrometer (A VIR1S) 

The next major step in the evolution of imaging spectrometers was the Airborne 
Visible/Inffared Imaging Spectrometer (AVIRIS). First proposed to NASA in 1983, 
AVIRIS first flew engineering evaluation flights in 1987 and became operational in 1989. 
AVIRIS was originally designed for the airborne collection of spectral images that would 
provide data necessary for the development of hyperspectral information extraction 
methods. The sensor is flown on a NASA ER-2 aircraft and can image in 224 contiguous 




spectral channels at a spectral resolution of 10 nanometers while covering the spectral 
region from 0.4 to 2.5 micrometers. Figure 2.8 shows the design layout of the AVIRIS 
sensor. 

Initial results from AVIRIS were promising. Image geometry, uniformity, 
signal-to noise ratio and spatial resolution were just as expected. The use of AVIRIS 
showed that valid spectral and radiometric calibration of radiance measuring instruments 
is required for physically based analysis of the measured data and for quantitative 
comparison of data acquired at different sites, and times. Figure 2.9 shows an early 
example of the spectral quality delivered by the AVIRIS sensor. This spectra was 
derived from averaging the individual spectra extracted from a 5 x 5 pixel region of an 
image acquired over the Cuprite Mining District in Nevada. It is interesting to note the 
agreement between the characteristic atmospheric absorption features that are observed 
with AVIRIS, (solid line) and the LOWTRAN atmospheric model predictions in 
parentheses. AVIRIS requirements initially specified a spectral calibration accuracy of 
+/- 5 nanometers. In-flight calibration of the spectral and radiometric characteristics of 
the sensor by Green et al, (1990) showed that a significant radiometric error could be 
introduced with as little as +/- 1 nanometer error in the wavelength calibration. The 
AVIRIS sensor has been continually upgraded since it's inception by improving the 
radiometric calibration and in-flight stability, increasing the signal-to-noise ratio, 
enhancing the scanner performance for geometric image fidelity and increasing the 
overall instrument reliability. 

Calibration experiments with the instrument have continued with each flight 
season. These periodic experiments combine AVIRIS overflights with in-situ 
measurements of the surface and atmospheric conditions present during the overflight. 
The measurements are used to constrain the radiative transfer modeling algorithm which 
predicts the upwelling radiance incident upon the sensor. The actual radiance measured 
by AVIRIS and the model predicted radiance are compared in order to validate the 
spectral and radiometric characteristics of AVIRIS in flight. Figure 2.10 shows a 
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comparison of the MODTRAN predicted radiance and the AVIRIS measured radiance 
acquired during a calibration/validation experiment over Ivanpah Playa, California on 7 
March 1991. Figure 2.11 shows the Ivanpah calibration/validation signal-to-noise 
performance of the AVIRIS along with the design requirements for the sensor. Table 2.1 
highlights the current AVIRIS characteristics along with past collection operations. 
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Numerous geological, botanical and atmospheric studies have been conducted 
using the AVIRIS instrument. The use of AVIRIS for the detection of changes in 
vegetation communities between seasons has shown dramatic spectral differences within 
various vegetative communities (Miller, et al, 1990). Figure 2.12 shows the representative 
reflectance curves attained for various vegetation communities on April 13, 1989, June 2, 
1989 and August 31, 1989. The plots show distinctive changes in the red-edge 
reflectance during the April to June period and again from the June to August period. 

Hyperspectral imagers have been used for the measurement of primary 
productivity in small lakes. This activity is commonly derived by measuring the quantity 
of photosynthesizing material (e.g. phytoplankton). AVIRIS images acquired over 
Mono Lake in 1989 provided reflectance spectra indicative of chlorophyll a, (Figure 
2.13). An attempt was made to calculate the amount of Chlorophyll contained in the lake 
by using a chlorophyll concentration algorithm developed for the Coastal Zone Color 
Scanner (CZCS), carried on the Nimbus series satellites. The amount of chlorophyll 
calculated was half that determined by on-site sampling. The discrepancy was 
determined to be in the algorithm itself since it was designed for much smaller 
concentrations normally found in oceanic waters. Further development is underway on 
an algorithm that can be applied to inland water bodies of this size (Melack and Pilorz, 
1990). 

Total column water vapor determinations can enhance studies in the fields of 
meteorology, climatology and high precision geodesy using the global positioning 
system. Gao and Goetz (1990), were very successful in their determinations of column 
water vapor estimations using AVIRIS data. By using spectral curve fitting of AVIRIS 
acquired spectra and simulated spectra along with band ratioing techniques, Figure 2.14. 
and Figure 2.15, they were able to estimate the column water vapor amounts to within 5% 
of those determined by a radiosonde. 


13 



c. Geophysical Environment Corporation Imaging Spectrometer (GERIS) 

One of the major commercial sources for spectrometers is GER. One of their 

primary systems is GERIS, a 63 channel spectrometer that utilizes three spectrometers. 
The optics scanning mirror scans an area of 512 pixels covering out to 45 degrees either 
side of the nadir crosstrack position. The incoming signal is then split and sent to the 
three spectrometers, each of which has it own line array. Spectrometer I records signals 
in 31 channels between 0.47 to 0.84 micrometers at a spectral resolution of 12.3 
nanometers. Spectrometer II records 4 channels in the 1.40 to 1.90 micrometer range at a 
resolution of 120 nanometers. The SWIR region, 2.0 to 2.45 micrometers, is covered by 
spectrometer 3 at a resolution of 16.2 nanometers or 28 channels. 

GERIS was utilized in a study of a com stand that overlaid a former waste 
deposit site in Germany. The study was part of the European Imaging Spectrometry 
Campaign (EISAC). The waste site had opened in 1899 and was a depository for 
industrial and household refuse until it's closure in 1949 when it was covered with humus 
and converted to agricultural land use. Figure 2.16 shows the average spectra from three 
different locations on the study site. In the visible and near-infrared regions, (top plot), 
differences in relative reflectance are 2% and 5% respectively. In the SWIR region the 
differences are smaller and discontinuous. Lehmann, et al, (1990), determined that the 
relative differences were most probably caused by heavy minerals which caused a stress 
situation to occur in the vegetation cover. 

d. Fluorescence Line Imager (FLI) 

The Canadian Department of Fisheries and Oceans operates the Florescence 
Line Imager (FLI). The instrument was built by Moniteq Ltd and can be operated in one 
of two separate modes. In the spatial mode the sensor operates as a pushbroom, 
multispectral scanner. In the spectral mode, the sensor operates as an imaging 
spectrometer, imaging in 288 contiguous bands in the region from 0.4 to 0.8 micrometers. 
Each spectral band is spaced at 1.3 nanometer intervals which provides a spectral 
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resolution of 2.6 nanometers. The sensor is normally flown on an Ontario Centre for 
Remote Sensing aircraft at an altitude of 1000 to 1800 meters. Rock, et al, (1990), 
utilized the FLI while studying coniferous forest damage on Whiteface Mountain, NY. 
On-site, red edge (660nm-800nm) spectral samples of birch leaves and Balsam Fir 
needles were acquired at the time of the experiment. First derivative curves were 
calculated and the mean, +/- the standard deviation was plotted as shown in Figure 2.17. 
The FLI was then flown over the site. After calibration of the data from instrument 
received radiance to ground reflectance values, both the reflectance and first derivative 
curves for the Birch and Fir spectra were plotted as shown in Figure 2.18 Comparison of 
the in-situ spectra and the FLI acquired spectra indicates a remarkable similarity in the 
red-edge region. 

e. Hyperspectral Digital Imagery Collection Experiment (HYDICE) 

HYDICE was initiated in 1991 as a classified study to determine the utility of 
using hyperspectral imaging technology for intelligence applications. The effort was 
expanded in 1992 to include the procurement of a research quality sensor system with 
civil agency participation under the Congressionally-funded, Dual Use Initiative. 

Hughes Danbury Optical Systems (HDOS) has produced a highly calibrated, 
imaging prism spectrometer of the nadir, push-broom type. HYDICE measures reflected 
solar energy along a 1 km ground swath at a flight altitude of 6000 meters. The 
instrument's spectral range is from 400 to 2500 nanometers, with a nominal spectral 
resolution of 10 nanometers. Ground sampling distance is dependent on sensor flight 
altitude and ranges from 1 to 4 meters. 

As of winter 1995 HYDICE was being fitted into a Convair-580 aircraft for 
initial aircraft compatibility checks prior to the initiation of the 1995 data collection 
experiments. Figure 2.19 is a physical view of the HYDICE instrument. Signal-to-noise 
ratio is an important requirement for earth remote sensing. Since the entire spectral range 
of the instrument is incident on the same focal plane array, HYDICE incorporates a 
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multiplexer to accommodate the large dynamic range required by the focal plane. As 
Figure 2.20 indicates, the dynamic range is broken down into three regions. A, B and C. 
The multiplexer is also divided into three regions to accommodate the maximum 
expected signal. This design effectively optimizes the signal-to-noise ratio together with 
the dynamic range. 

A spacebome derivative of the HYDICE concept has been investigated. 
(Silvergate, et al, 1994). By using telescope type foreoptics the HYDICE concept could 
be adapted to a spacebome configuration, see Figure 2.21. 

f Shuttle Imaging Spectrometer Experiment (SISEX) 

Designed as a demonstrator of technologies relating to imaging spectrometers, 
SISEX would image in 128 separate spectral channels from 0.4 to 2.5 micrometers, 
provide a ground instantaneous field of view (GIFOV) of 30 m at a Shuttle altitude of 
approximately 250 km, and produce a swath width of 12 km. Spectral resolution in the 
early stages of the design was set at 20 nanometers but was sharpened to 10 nanometers 
after field work and analysis in support of the AVIRIS concept determined the latter as 
more beneficial (Herring, 1987) After the design process began, the Earth Observation 
System (EOS) was conceived. Initial EOS plans called for a High Resolution Imaging 
Spectrometer (HIRIS) to be deployed with the system. In an effort to conserve dollars, 
NASA began to develop SISEX as a test bed for the HIRIS design and planned to utilize 
SISEX components on HIRIS in order to reduce costs. The Challenger accident caused 
NASA to restrict any future use of the shuttle for advanced instrumentation testing and in 
1987 the Jet Propulsion Laboratory canceled SISEX in favor of HIRIS (Taranik, et al, 
1993). 

HIRIS development continued and in 1988 the Jet Propulsion Laboratory 
proposed to fly HIRIS on a polar orbiting platform. HIRIS was to image in 196 spectral 
channels in the range from 0.4 to 2.5 micrometers. However, the declining NASA budget 
caused an indefinite postponement of the HIRIS program in 1992. 
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g. Moderate Resolution Imaging Spectrometer (MODIS) 

MODIS is an Earth Observing System (EOS) instrument that is slated for flight 
on board the EOS-AM and EOS-PM satellite series. MODIS will be employed to enable 
a better understanding of the global earth system to include the interaction between land, 
ocean and atmospheric processes (Thompson, 1990) The instrument is an imaging 
radiometer incorporating a traditional, cross-track, scan mirror, coupled with linear 
detector arrays. There are four focal planes that contain the detector arrays and each 
incorporates spectral interference filters. 

MODIS will image in 36 spectral bands in the range from 0.4 to 15 micrometers. 
Spatial resolution will be between 250 meters and 1 kilometer at nadir. 
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3. Evolution Summary 

The following table summarizes the evolution of hyperspectral imaging sensors. 
The Landsat Thematic Mapper is included to enable the reader to compare the attributes 
of current multispectral sensors with those of the hyperspectral breed: 



Table 2.2 ( Bandwidth characteristics of multispectral and hyperspectral systems) 


B. ANALYSIS TECHNIQUES 

Multispectral image analysis techniques can further exploit information contained 
in multispectral images. Visual classification can be achieved by assigning three of the 
primary colors (red, green, blue) to each of the three sensor bands selected and viewing 
the image. The information content represented by this false-color image will depend on 
the spectral bandpasses utilized. 
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Land cover categorization, more commonly known as Terrain Categorization 
(TERCAT) involves categorizing the image pixels into land coverages utilizing spectral 
pattern recognition. Figure 2.22 shows a 3-dimensional space defined by the axes which 
represent the intensities of the selected bands, in this case, 3 bands from the LANDSAT 
Multispectral Scanner (MSS). A point can be defined within this space for each pixel in 
the image. In an ideal case, each class or category of land cover would be represented by 
the same point. In reality, the category will be represented by a clustering of pixel values. 
The different clusters may be distinct and widely separated from each other and can be 
identified with a specific land cover type. This generally requires that a study of the area 
in question has been done to determine the predominant land cover types present. 
Otherwise, simple operator experience can be used as a basis. By using selected pixels 
from the scene, commonly referred to as "training data" or "sites", high levels of 
classification accuracy can be achieved. A number of techniques have evolved which 
constitute rotations in wavelength (color) space. 

1. Principle Components 

Principle components or the Karhunen-Loeve transformation grew out of the 
statistical treatment of psychological problems (Cracknell and Hayes, 1991). Principle 
components is a special data transformation which transforms the information content of 
the original measurement set (e.g., multispectral bands) into a coordinate system which is 
orthogonal in it's new coordinates. This orthogonality is determined by a rotation that 
leaves the data uncorrelated in the new coordinates. In some instances, these new 
principle component bands may be more easily categorized than the original image. 

Principle components analysis has the following advantages in the analyses of 
multispectral and hyperspectral data: 

1) Most of the radiance in the multispectral data is relegated to one or two 
principle component bands. 

2) Noise relegated to less-correlated principle component bands. 
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3) Spectral differences between different surface materials become more 
apparent in the principle component bands than in the original multispectral 
bands. 

Application of principle component transforms to hyperspectral data has only 
recently become practical. It will help in categorization, and reduction of the 
dimensionality of the data sets. The number of principle component images generated 
will be equal to the number of spectral bands that make up the image cube. The majority 
of the information within the scene will most likely be described by the first few 
(generally 10 or less) principle components. Principle components will reduce the data 
volume of the hyperspectral set but will de-emphasize individual spectra or target 
signatures. This is due to the fact that the majority of the spectral signatures within a 
hyperspectral scene are highly correlated not orthogonal, and the principle component 
images contain information that is a linear combination of several spectral classes 
(Harsanyi and Chang, 1994). 

Sub-pixel targets can only be identified by utilizing data reduction techniques that 
model the spectral contributions of the target and background materials (Farrand and 
Harsanyi, 1995). Spectral mixture analysis techniques and low probability of detection 
algorithms based on orthogonal subspace projection show the greatest promise for 
detection and classification. 

2. Spectral Mixture Analysis 

Spectral mixture analysis is based on the assumption that the majority of the 
spectral variation in a multi or hyperspectral data set can be accounted for by a relatively 
small number of endmember spectra. Endmembers define the vertices of the mixing 
space. The first step in the spectral mixture analysis process is to determine the number 
and identity of the spectral endmembers. Image analysis software such as the Spectral 
Image Processing System (SIPS), developed by the Center for the Study of Earth from 
Space (CSES), utilize programs written in Interactive Data Language (IDL) format for 


20 



linear spectral unmixing. During the process, a spectral endmember library is formed 
using an iterative process and then decomposed using Singular Value Decomposition 
(SVD). This allows for determination of the degeneracy of the library or the number of 
eigenvectors that make the largest contribution to the variance spanned by all the 
endmembers. One method of determination has been the noise-adjusted, Principle 
Components Transform. If the results of the degeneracy indicate that the normalized 
singular values of the endmembers are unique (spectrally separable and orthogonal), then 
the user can continue on to the next phase of the unmixing process. However, if all but 
one of the singular values is zero, the user must reselect new endmembers prior to 
proceeding. The endmember spectra can then be extracted directly from the Principle 
Component images and will define areas on the ground that are spectrally unique and 
which correspond to specific materials. Laboatory or field measured spectral 
reflectances can then be used in the spectral mixture analysis program as "reference 
endmembers". The user can then run the unmixing program whose output is another 
image data cube. The cube contains one image for each endmember showing the derived 
spatial patterns of abundance for that endmember. Additionally, two additional images 
are generated that can be used to assess the uncertainty in the unmixing results. 

3. Low Probability of Detection (LPD) 

LPD is based on the concept of orthogonal subspace projection (OSP) which is a 
result from least squares theory (Harsanyi and Chang, 1994). The technique basically 
involves a two step process. The first step is the determination of a matrix operator that 
eliminates undesired or interfering signatures. Basically, this is an optimal interference 
rejection process in the least squares sense. Secondly, it is necessary to develop a vector 
operator that maximizes the residual desired signature signal-to-noise ratio. Both of these 
operators are then combined into an overall orthogonal subspace projection classification 
operator which reduces the non-Gaussian detection and classification problem presented 
by pixels containing mixed spectra to the solved problem of detecting an unknown 
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constant in white noise. The LPD algorithm does not require the user to have prior 
knowledge of the background signatures present in the image. In a simulation conducted 
by Harsanyi and Chang, 100 mixed pixels were simulated using reflectance spectra for 
red soil and dry grass and creosote leaves, as shown in Figure 2.23. Pixels 20, 40, 60 and 
80 contained the creosote leaves reflectance spectrum at abundances of 20, 15, 10 and 5% 
respectively. Figure 2.24, illustrates the detection and classification problem associated 
with spectrum of similar shape. In this case, pixels 18-22 are shown on the plot. Pixel 
20, which contained the creosote leaves reflectance spectrum at 20% abundance is 
included in the plot but is not readily distinguishable from the other pixels. The OSP 
operator was then formed using the creosote leaves spectrum as the reflectance signature 
of interest. The operator was then applied to each pixel. Figure 2.25 indicates that the 
pixels containing the creosote leaves spectrum were selected by the algorithm. 

The application of this technique to hyperspectral images with mixed pixels has 
shown that signatures of interest can be detected at abundance levels as low as a few 
percent at signal-to-noise ratios of less than or equal to 50:1 (Harsanyi and Chang, 1994). 

4. Constrained Energy Minimization (CEM) 

CEM is an additional analysis technique that is similar to LPD. However, the CEM 
operator does not depend on the signature of interest occurring with a low probability. It 
is less stringent than the LPD approach identified above and like LPD, does not require 
prior knowledge of the undesired background signatures. In basic terms, a linear operator 
is used that minimizes the total energy in a hyperspectral image sequence with the 
constraint that the energy associated with the signature of interest is maintained. In this 
case, CEM works well in the detection of sub-pixel targets that occur over a large number 
of pixels. The CEM technique is particularly useful for geologic and environmental 
research as demonstrated by Harsanyi et al, (1994). 
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III. EXPERIMENT DESCRIPTION 


The data utilized in this thesis is the result of airborne hyperspectral remote sensing 
collections conducted at White Sands, New Mexico in October 1994. The name of the 
experiment was Operation DESERT RADIANCE. 

The HYMSMO program sponsored a comprehensive survey of military 
operational and intelligence needs. The needs are currently under evaluation by the 
Unified Commands for prioritization based on time-sensitivity and information criticality. 
In addition, the Defense Intelligence Agency (DIA), sponsored a study that determined 
the information needs of military commanders based on operational user definitions. 
These "levels of information" (LOI’s) show common areas that can be used to effectively 
reduce the HYMSMO needs to a common set of experiments. These experiments are 
based on target classes recognized by the Joint Chiefs of Staff (JCS): objects (mobile and 
fixed), emissions, activities and events. 

A. OPERATION DESERT RADIANCE 

The DIA sponsored, Assured Support to Operational Commanders (ASOC) study, 
identified LOI's which include: characterization (target spectral "signatures"), detection 
classification, discrimination and identification (material, function and unit). Each of 
these LOI's relate to the set of target classes previously identified. 

The DESERT RADIANCE experiment sought to concentrate it's efforts in the 
LOI's concerned with exposed objects (mobile and fixed) and cultural features. The main 
reasoning for this was the necessity to gain experience and knowledge from exposed 
object signatures in order to properly evaluate the degradation of LOI content as the 
target becomes obscured or camouflaged. The underlying principle being that the ability 
to characterize the spatial, spectral (textural, temporal and other) signature of a target is 
necessary for detection and satisfaction of the LOI's. DESERT RADIANCE'S highest 
priority was to determine the ability to remotely detect and identify the material 
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properties of a target at the subliteral and/or subpixel level, given proper spectral 
characterization. Abundance estimation or the quantification of a target material type at 
the subpixel level and evaluating the effect of spectral taggants were additional 
objectives. To further define the target class of exposed mobile and fixed objects 
utilized for this experiment, one must consider that mobile objects are generally those 
that are less than or equal to a few pixels in the spatial dimension which would be 
adequate for a literal imagery analyst to perform detection and geolocation based on 
spatial characteristics alone. The fixed objects and/or cultural features are those that 
occupy more than a few pixels and are termed "multi-pixel" targets. 

1. Experiment/Objective Breakdown 

The DESERT RADIANCE operation involved objectives relating to the airborne 
collection of target signatures and objectives relating to the tower-based collection of 
target signatures. 

a. Airborne Data Collection Objectives 

Two experiments involved airborne data collection. Experiment A dealt with 
exposed mobile objects located at test sites A and B. Experiment B examined exposed 
fixed objects, cultural features located at White Sands and Las Cruces, N.M. The 
following is a list, by experiment, of the objectives: 

(1) Experiment A 

♦ Remote spectral characterization of materials associated with exposed 
mobile objects. 

♦ Remote subliteral/subpixel detection of separate categories of materials 
associated with exposed mobile objects. 

♦ Remote identification of materials associated with exposed mobile 
objects, at both the literal and subliteral level of spatial resolution. 
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♦ Remote detection and functional identification of exposed mobile 
objects, at both the literal and subliteral level of spatial resolution, through 
the use of spectral taggants. 

(2) Experiment B 

♦ Remote spectral characterization of materials associated with exposed 
fixed objects. 

♦ Remote identification of materials associated with exposed fixed 
objects, at both the literal and subliteral level of spatial resolution. 

♦ Remote detection, classification and discrimination of separate 
categories of materials used in exposed fixed objects and/ or cultural 
features. 

♦ Quantification of the aerial extent and relative abundances (subpixel 
fraction) of materials used in exposed fixed objects and/or cultural 
features. 

♦ Remote detection and functional identification of exposed fixed objects, 
at both the literal and subliteral level of spatial resolution, through the use 
of spectral taggants. 

b. Tower-Based Data Collect Objectives 

The tower-based data collect concentrated on the spectral characterization of 
exposed mobile objects and desert backgrounds. The objectives are as follows: 

I) Assess the variation in spectral signatures as a function of solar angle: 

* Complex target, full pixel 

* Simple target, partial pixel 

* Complex target, partial pixel 
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2) Remote spectral characterization of desert backgrounds. 

3) Remote spectral characterization of exposed mobile objects at long 
range. 

2. White Sands Missile Range (WSMR) 

The WSMR is located 40 miles north of El Paso, Texas in southern New Mexico. 
Figure 3.1 depicts the relative location of the WSMR. The main test area is 
approximately 40 by 100 NM in size and is generally utilized by the Army to evaluate 
missile and rocket systems. 

The primary site of interest at WSMR for the operation was the Electro Optical 
(EO) Range. Located within the range is a 100' tower, experiment test site A and 
experiment test site B. The two experiment test sites were developed to take advantage of 
different vegetative backgrounds associated with each site. Each test site was combined 
with the 100' tower location to form separate flight tracks followed by the airborne 
sensors. In addition, Condron Airfield, located to the south of the tower was used as an 
additional calibration/visual checkpoint for experiment site A. Figure 3.2 depicts the 
actual flight patterns flown for the data collections on 27 and 28 October 1994. 

3. Target Descriptions 

a. Fabric and Painted Panels 

The Topographic Engineering Center (TEC) supplied the fabric and painted 
materials to facilitate the construction of target panels. These panels were grouped in 
some cases to form target panel arrays that simulated exposed mobile object dimensions. 
Panel dimensions were varied in size to allow some panels to either completely fill pixels 
or be subpixel depending on sensor altitude. These dimensions are based on the expected 
ground sample distance (GSD) calculations as well as the along-track and lateral 
instantaneous fields of view of the sensors. In addition, panels of similar material were 
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placed close together to allow for spectral mixing when viewed from high altitudes. 
Extreme care was taken with the target panels. Each night, all panels were removed from 
the experiment test sites and covered to prevent any moisture/dust contamination when 
not in use. WSMR personnel supplied two camouflaged nets that were used to conceal 
vehicles and desert background materials, fable 3.1 lists the fabric materials and Table 
3.2 lists the painted panels that were deployed during the DESERT RADIANCE 
exercise. 



Table 3.1 (Fabric Materials Deployed, Operation Desert Radiance, Anderson,1995a) 
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Designation/ 

Priority 

Target Target 

Description Size (ft) 

Mix 8 '* 

Pattern 

pi /1 

Black 

10x6 

c4 

P2 /I 

Brown 

10x6 

c4 

P3 n 

Tan 

10x6 

c5 

P4 / 1 

Tan Taggant (SpectIR) 

10x6 

c5 

P5 /I 

Green 

10x6 

c4 

P9 /I 

Lt Tan Taggant Surface Optics) 

10x6 

c5 

P10/1_ 

Lt Tan Reference(Surface Optics) 

10x6 

c5 


Table 3.2 (Painted Panel Deployed, Operation Desert Radiance, Anderson, 1995a) 


The WSMR provided a limited number of military vehicles and equipment to 
support the experiments. Table 3.3 lists the vehicle designations and descriptions. 


Designation / 
Priority 

Target 

Description 

Target 

Size (ft) 

!%-, 

VI /I 

camouflage 

9x20x8 

With/without F9 

V2/1 

M-60 tank, woodland 
camouflage 

9 x 20 x 8 

WilWwithout F10 

V5 /I 

Bronco, woodland 

5x8x4 

Open and near bushes 

V6/1 

Bronco, woodland 

5x8x4 

Open and near bushes 

V7 /I 

Bronco, woodland 

5x8x4 

Open and near bushes 

V8 /I 

Bronco, woodland 

5x8x4 

Open and near bushes 

V9 /2 

HMMWW decoy 

7x15x6 

Open and near bushes 

V10/1 

CHAPARRAL 

8 x 15x8 

In tower experiment 


Table 3.3 (Tactical Vehicles Deployed, Operation Desert Radiance, Anderson, 1995a) 
Calibration targets were deployed by EG&G, Incorporated. The targets were 
used to provide reflectance and absolute radiance references for the airborne data. Table 
3.4 lists the calibration targets and their associated reflectance percentages. Note that the 
60% reflectance target was not deployed due to the likelihood that it would saturate the 
AAH1S airborne sensor. Figure 3.3 depicts the Site B target layout for 27 October 1994. 
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b. Optical Target 

An optical test target was deployed during the collection runs. It's design was 
based on a concept developed by Dr. William Stoner, Science Applications International 
Corporation (SAIC). Dr. Stoner proposed that the exposure of small reflective gratings 
to sunlight will produce diffracted light at discrete wavelengths depending on the 
orientation of the gratings (Stoner, 1994). 

Figure 3.4 depicts the geometry of the sun, grating and airborne sensor. The 
grating grooves are perpendicular to the plane of the sun, grating and sensor flight path. 
This configuration presents the possibility that discrete line spectra can be produced and 
detected by the sensor. Figure 3.5 shows the grating rotated from the perpendicular to 
the sun, grating and sensor flight path. By comparing the two figures one can see that the 
reflected light rays are unaffected by the rotation, but the diffracted rays are rotated 
within the plane of the sun equal to the rotation of the grating. The diffracted ray that is 
"up and normal" to the grating would be detectable by the airborne sensor flying within 
the same plane. Rotation of the grating will cause a cone to be formed by the diffracted 
rays about the reflected ray, see Figure 3.6. 

Dr. Stoner proposed that circular gratings in the form of compact discs (CD's) 
would provide the necessary grating geometry to produce a nearly continuous cone of 
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diffracted rays about the reflected rays and accommodate uncertainties in the knowledge 
of the sun's position. 

The solar spectrum is fairly smooth at the 10 nanometer scale. The collected 
wavelength spikes are however more narrowband in nature. This is due to the fact that 
the geometry of the sun, grating and sensor functions like a grating monochrometer. The 
frequency of the wavelength spikes will be dependent on the sun's zenith angle and 
azimuth which can be derived for specific locations given the time of day. 

The total spectrum collected by the sensor will be a combination of three 
components: 

1) Line spectrum from deffracted sunlight. 

2) Skylight reflected specularly by the gratings. 

3) Light reflected from the natural background in the remaining portion 
of the ground-spot that is not taken up by the gratings. 

Relative radiances were calculated to determine an optimum number of CD's for 
the experiment. Four CD's were utilized and placed on a self-leveling platform elevated 
10 degrees. Each CD was then placed on the platform at 0, 2, 4 and 6 degrees elevation 
respectively to provide varying orientations with respect to the sun. This set-up was 
necessary to provide four separate wavelength spikes during the sensor overflight. 
Wavelength spikes can be determined by using a Mathcad program developed by Dr. 
Stoner and observations of sun zenith angle and elevation taken at the site starting at 
1000 (L) through 1600 (L), 27 October 1994. 

4. Ground Truth Measurements 

The Topographic Engineering Center (TEC) was responsible for formulating the 
Ground Truth Protocol for the operation. The Protocol defined the standard operating 
procedures for collecting ground truth measurements in the field. The purpose of the 
ground truth collection during the airborne collection operations was to provide a means 
to relate the reflectance of the surface to the energy recorded by the imaging sensor. The 
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energy reaching the sensor is dependent on the incident irradiance at the surface, the 
physical properties of the surface, emission of the surface and the effects of atmospheric 
absorption and scattering on the incident and reflected energy. In addition to TEC's in 
situ spectroradiometric measurements of the target materials and background, EG&G 
provided ground truth measurements of the calibration target panels. The spectral data 
was compiled and made available as an ascii file to facilitate any spectral comparison of 
the target/background/calibration materials with image derived spectra. 

5. Sensor Characteristics 

The DESERT RADIANCE experiment set forth specific requirements for sensor 
performance. Availability and cost were the initial drivers for the selection process. 
However, due to the nature of the experiment (detection/descrimination of targets from 
their background), and it's hyperspectral orientation, several key sensor performance 
characteristics needed to be met. The following characteristics were considered in the 
final selection: signal-to-noise ratio, number and width of the spectral channels, overall 
spectral coverage range and adequate spatial resolution. As a result of a technical 
assessment of the available sensors, two sensors were selected for the DESERT 
RADIANCE experiment. The SETS Technology, Advanced Airborne Hyperspectral 
Imaging System (AAHIS) and the Geophysical and Environmental Research (GER) 
Corporation's, 63 channel, Digital Airborne Imaging Spectrometer (DAIS). 

a. Advanced Airborne Hyperspectral Imaging System (AAHIS) 

The AAHIS sensor is a second generation sensor that was originally developed 
by Science Applications International Corporation (SAIC), San Diego, California. The 
sensor was upgraded in a cooperative effort with SETS Technology INC., Honolulu, 
Hawaii and is utilized by the state of Hawaii for maritime applications. The basic block 
diagram depicting the AAHIS major components is shown in Figure 3.7. The AAHIS 
sensor was designed for high signal-to-noise ratio performance based on it's intended use 
in maritime applications which involve low reflectance environments (e.g. 5% 
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reflectance). Table 3.5 lists the AAHIS performance characteristics. The table also lists 
the sensor's signal-to-noise ratio performance over a land-based (20% reflectance) 
environment. Figure 3.8 shows the signal-to-noise ratio performance in both 
environments over the full wavelength span of the sensor. The sensor was calibrated at 
the Naval Command, Control and Ocean Surveillance Center, Research and Development 
Division (NRaD), using a 1 meter integrating sphere and subsequently installed into a 
Piper Aztec, twin engine, aircraft for flight operations at WSMR. 


Sensor Perimeters 

Nominal Value 

SignaJ-to-Noise (5 % reflectance target) 

■ 50 Hz frame rate 

200-400 

Signal-to-Noise (20 % target) 

• 50 Hz frame rate 

500-850 

Useful spectral range (nm) 

440-870 

Number of spectral channels 

72 

Nominal bandwidth (nm) 

• 4 channels summed on chip 

• 2 channels summed in software 

12.4 

Cress track IFOV (mradians) 

i.t 

Along track IFOV (mradians) 

1.0 

Spatial pixel number 

190 

Frame rate (maximum) (Hz) 

55 

Swath width (degrees) 

11.4 

Data record length/tape at 55 Hz (seconds) 

30 

Digitizer (Bits) 

12 


Table 3.5 (AAHIS Sensor Characteristics, Anderson, 1995a) 


b. GER Digital Airborne Imaging Spectrometer (DAIS) 

The GER sensor was selected to compensate for the lack of a short-wave 
infrared (SWIR) capability in the AAHIS sensor. Many unique spectral features manifest 
themselves in the SWIR domain in narrow spectral bands thereby validating the need for 
a sensor such as the GER DAIS. Figure 3.9 depicts the sensor configuration and optics 
layout. The DAIS employs 3 different detector arrays, Silicon for the VIS/NIR, Lead 
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Sulfide (PbS) for the SWIR and Mercury Cadmium Telluride (HgCdTe) for the thermal 
IR The sensor's low scan rates can affect the spatial resolution at low altitudes and cause 
subsequent undersampling of the scenes in this regime. However, the sensor can still 
provide high spectral resolution in the same context. fable 3.6 lists the sensor 
characteristics. Calibration of the DAIS sensor is not possible due to the wide scanner 
aperture and the lack of a calibration facility, therefore a 50' x 50’ area of the Condron 
Airfield runway approach end was used as a calibration site. The selection of this site 
was based on its relative uniformity with respect to background. Ground truth data of the 
calibration site was taken by TEC. 



Table 3.6 (GER 63 DAIS Sensor Characteristics, Anderson, 1995a) 
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IV. EXPERIMENT OBSERVATIONS 


Observations of the experiment conducted by the HYDICE office and exploitation 
team is presented below. This chapter is broken down into twomain sections. The first 
presents the observations made by the author of the experiment test sites, target 
preparation and handling, and sensor platform performance during the set-up and conduct 
of Operation Desert Radiance at White Sands Missile Range (WSMR) from 24-25 
October 1994. The second section is the observations of the hyperspectral data sets 
produced by each of the collection platforms. 

A. WHITE SANDS MISSILE RANGE OBSERVATIONS 

The White Sands Missile Range was selected for the experiment based on the 
spectral consistency of the desert background present, the predictability of weather 
conditions along with low humidity, tower access and general capability of the Missile 
range complex to accommodate and facilitate an exercise of this type. The facility is 
located at an altitude of approximately 4000' and lies in a relatively flat valley on the 
eastern side of the Organ Mountains in southeastern New Mexico. The valley is oriented 
N-NE to S-SW with a gradual elevation increase as one travels northeast. 

1. OPERATION DESERT RADIANCE 

The DESERT RADIANCE experiments were to focus on overtly exposed mobile 
and fixed targets, to include cultural features and backgrounds. The author observed 
experiment set-up and participated in the preparation of targets and target placement sites 
on 24-25 October 1994. The actual collection flights occurred after the author's departure, 
on 26 , 27 and 28 October 1994. Therefore, the author does not have first-hand 
observations relevant to these collection periods. Minor adjustments to portions of the 
experiment set-up occurred during these two periods, but they do not appear to have 
affected the outcome of the experiments. 
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The experiment effort was planned to take advantage of the exploitation team's 
manpower. Separate groups were formed within the exploitation team to facilitate target 
construction, experiment site preparation and other specific technical duties such as 
ground truth data collections. 

2. Site Layout 

The White Sands Missile Range is uniquely suited to accommodate an experiment 
of this type. It's isolation from major urban areas and overlying restricted airspace 
provide a controlled environment within which airborne sensor platforms can maneuver 
and the desert environment provides a spectrally uniform environment that can enhance 
the contrast between target and background. 

Within the range complex is the Electro Optical (EO) Range that contains a 100' 
tower, clear dirt pad areas and well maintained, dirt/gravel access roads. The composition 
of the roads enhances the production of dust clouds whenever a vehicle passes but the 
two experiment sites appeared to be at a sufficient distance from this contamination 
source. 

The two experiment sites, A and B, are within the EO Range and can be observed 
easily from the tower. Site A appears to run N-NE from the tower and tower access road 
on the southern end to an access road running NW-SE on the northern end. Site B is 
approximately 1/4 mile NW of the tower and contains a 50' wide (approx.) dry wash and 
is bounded by an access road on the NE side. 

The tower itself has large landings every twelve feet with 360 degree look 
capability. The top of the tower incorporates a covered structure which contains 
electrical outlets, tables and a small, accessible balcony for equipment mounting. 

3. Weather Conditions 

The southeastern portion of New Mexico provides a benign weather environment in 
which to conduct experiments. The latitude of the site provides for a relatively high sun 
angle, with very little influence from frontal activity year round. The overhead 
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conditions on 24 October 1994 ranged from clear skies to a high percentage cloud 
coverage, with altocumulus and altocirrus cloud types present. The first flight day, 25 
October, held marginal weather conditions consisting of cumulus clouds with bottoms at 
5000' above ground level (AGL). Temperatures varied from the 40's overnight to the mid 
70's during the daytime. Relative humidity was low, but the low overnight temperatures 
caused dew to form on the landscape. Wind was usually calm in the morning hours but 
began to pick up slightly in the early afternoon. 

4. Targets 

The targets utilized in OPERATION DESERT RADIANCE were detailed in 
Chapter III. The WSMR provided workshop facilities for the construction of target 
frames and the mounting of target materials. Upon arrival at the target construction site, 
the author observed numerous targets in various stages of completion. Distinct target 
sizes were observed and closely approximated the target dimensions identified in the 
previous chapter. Tactical vehicles with various paint schemes were utilized. As these 
vehicles were moved along the dirt/gravel roads, dust was thrown up. No apparent 
attempt was made to clean the dust collected on the vehicles prior to any overhead 
collect. 

a. Materials 

Materials for the targets were provided by the Topographic Engineering Center 
(TEC) and brought to the experiment site. The fabrics appeared to be those used in 
military clothing, tenting, parachutes and camouflage, (e.g. Nomex, desert camouflage 
material, woodland camouflage material and green nylon parachute material). The color 
make-up of the fabrics included dark/light greens, tans, tan/browns and black. The 
parachutes were not homogeneous in color but instead, were made up of triangular 
segments with varying shades of Olive Drab (OD) green. The common thread was that 
all the fabrics were designed to inhibit detection by blending in with surrounding 
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background materials. TEC also supplied 2' x 2' Aluminum panels painted with various 
types and colors of paint. 

b. Construction 

Each target panel was constructed with 1" x 4" Pine boards with each 
connection secured with wood screws. A hole was drilled into each comer of the panel 
frame to accommodate 10-12" spikes used for anchoring purposes. The individual fabric 
types were then cut to dimension and stapled onto the frames. The size of the frames and 
the composition of the fabrics limited the amount that the material could be stretched 
across the frames without tearing. This caused the fabric to dip slightly as it layed across 
the frame. In addition, some warping of the wood frames was observed. Similar frames 
were constructed for the painted panels which were placed inside the frame structures to 
form 6' x 10’ target arrays. The individual panel comers were not secured and some of 
the panels were observed to warp while laying inside the frames. 

Parachute material was used to cover a 4’ x 4' x 4' box target and was folded 
to form various sized panels (without frames) and staked to the ground. The color of the 
parachute material stretched over the box was not homogeneous and was fairly 
transparent. A similar condition occurred with the parachutes folded and staked to the 
ground. The edges of the panels, which contained most of the folded material, appeared 
darker than the material towards the center of the panel. The parachute material could be 
stretched but was susceptible to tearing at the staked comers. In addition, the thinness of 
the material allowed for a small amount of material flapping during windy conditions. 

c. Care and Maintenance 

The exploitation team was provided with a large flatbed tractor trailer to 
facilitate the transportation of the target frames to and from the experiment sites. Care 
was taken in the loading/off-loading of the targets to eliminate tearing and minimize 
dust/dirt contamination. Due to the probability of moisture formation on the targets 
during the night and the possibility of particle contamination during transport, large 
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plastic sheets were secured over the top of the targets after they were stacked aboard the 
flatbed. 

d. Deployment/Retrieval 

Prior to each day's collection event, the target panels were deployed to the 
designated experiment site, A or B. Exploitation team members were formed into 
working groups to deploy/retrieve the targets from the sites. WSMR personnel were 
utilized in the movement of tactical vehicles to designated locations and the erection of 
camouflage netting. 

5. Experiment Sites 

Two experiment sites were utilized for the DESERT RADIANCE operation. Both 
sites were in close proximity to the observation tower and initially looked no different 
from one another. Closer inspection reveals that the sites provide separate background 
types. 

a. General Soil and Vegetation Make-up 

Site A contains a large number of Mesquite and Creosote shrubs ranging in size 
from 1 to 5 feet. The soil conditions are a loose sand and dirt mixture with small areas of 
hard-packed soil. The large quantity of shrubs enable targets to be placed adjacent to the 
shrubs thereby providing some coverage and spectral mixing. Site B contains much less 
vegetation and a dry wash on the southern end. Very loose sand makes up the wash and 
it's adjacent banks while more hard-packed dirt and sand encompass the remaining 
portion of the site. 

b. Preparation 

Extensive preparation of the two experiment sites was performed prior to target 
placement. Each placement location was marked out to the approximate dimensions of 
the target to be placed there. The area within the target dimensions was cleared of all 
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low-lying shrubs, grasses and debris and the soil leveled to ensure that once the target 
was placed, it presented a relatively level surface to the sensors. Care was taken to 
ensure that a minimal amount of soil was disturbed beyond the boundaries of the targets. 

6. Ground Truth Measurements 

Ground truth measurements were taken on a variety of objects both natural and 
man-made. Prior to the deployment of man-made target materials, ground truth 
measurements were collected on the surrounding desert vegetation and soils. In addition, 
the TEC teams used cherry pickers to collect spectra from the M-60 tanks and other 
tactical vehicles. The man-made material spectra were acquired by positioning the tripod 
of the portable field spectrometer adjacent to a target panel. The spectrometer was swung 
over the target material which positioned the instrument aperture approximately 10-12" 
from the edge of the target frame. The composition of the different fabric materials 
restricted the amount of stretch that could be applied without tearing. This fact caused, in 
certain instances, the material within the field of view of the spectrometer not to be 
perpendicular to the sensor. Observations of the calibration panels staked out by EG&G 
from the tower showed significant deformations in the panel surfaces. The ground truth 
data was recorded on portable PC’s from the spectrometers and copied in a DOS format to 
3.5" floppy disks. The number of spectra taken in the 0.4-2.5 micron range varied for 
each target of interest. 

7. Sensor Platform Performance 

The author observed one day of sensor overflights on 25 October, 1994. A flight 
planning session was conducted separately each night after the overall planning meeting 
to address issues or procedural problems that had arisen during the previous days flights. 
Collection altitudes were discussed and chosen based on the frame rates of the sensors 
and the desired GIFOV. It appeared that the flights at the higher altitudes of 4, 6 and 
8000' (AGL) would be affected by the pilot's inability to reference the flightline markers 


40 



layed out at either end of the track. The GER pilot visited the experiment sites and used a 
hand-held Global Positioning System (GPS) receiver to record navigational waypoints. 

a. Pilot Proficiency 

It became apparent during the first meeting on 23 October, 1994, that both pilots 
were relatively inexperienced. The author has no knowledge of the number of flight 
hours possessed by each pilot or their experience level in relation to remote sensing 
overflights. The pilots were concerned about the restricted area course rules, 
check-in/check-out procedures, experiment flightline orientations, flight altitudes, high 
altitude aircraft performance and bad weather flight determinations (weather calls). 
Concerns about the presence of two aircraft within the operating airspace of the 
experiments at the same time were alleviated after discussions with the experiment 
coordinators. 

Both aircraft attempted to fly the pre-briefed flight tracks on 25 October 1995 but 
they did not orient themselves properly. Both aircraft ran patterns from north to south 
instead of from south to north as briefed. The effect of crosswinds was very evident from 
the ground perspective. The ground track of the sensor platforms progressively moved 
east/southeast within the experiment area boundaries. This caused lateral displacement 
from the run-in course line, poor sensor-target geometry or a complete miss of the target 


b. Airspace Coordination 

Prior coordination with the WSMR range control was made during advance trips 
and prior to the actual flights. The initial check-in and entry into the restricted area did 
not go smoothly and some course rules were not briefed to the experiment team by 
WSMR personnel. Radio problems with one of the aircraft, frequency assignments and 
proper execution of course rules were some of the problem areas. No face-to-face briefs 
were conducted between range personnel and the sensor platform pilots. 
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c. Sensor Reliability 

The AAHIS sensor developed inverter problems during the first overflight on 
25 October, 1994. The aircraft returned to the Las Cruces airport and the part was 
ordered and replaced prior to the 26 October collections over Las Cruces. Saturation of 
the detector array at high altitudes was a concern of the experiment team. 

There had been concern expressed regarding the frame rate of the DAIS sensor 
prior to the operation. GER stated that the sensor would provide a frame rate of 12 Hz. 
During the first coordination meeting on the evening of 23 October, 1994, the GER team 
stated that the sensor could only provide a frame rate of 8.5 Hz. This caused a concern 
over the fact that scenes imaged at the lower altitudes would be undersampled due to the 
aircraft speed over the ground. 

The GER aircraft experienced mechanical difficulties after the 25 October flight 
and was not utilized again until 28 October, 1994. 

B. HYPERSPECTRAL OBSERVATIONS 

Data from both sensors was provided to the author by the HYDICE Office through 
Dr. Linda Kalman, The Aerospace Corporation. In all, seven Exabyte tapes were 
received. Due to the limited disc space and the large volume of data contained on the 
tapes, a decision was made to load only certain files that had been found by the 
exploitation team to be of sufficient quality to be used in analysis work. An initial 
viewing of the AAHIS collected, 27 October data 8 file was conducted during late 
January and February 1995. Subsequent analysis was conducted at the HYDICE 
Program Office. The HYDICE exploitation team helped/directed the analysis, and as a 
result, the time to analyze the data was reduced significantly. 

1. Hardware and Software Configuration 

The hardware configuration at the HYDICE office consisted of a UNIX based, Sun 
workstation with approximately 4 Gigabytes of storage capacity. The software utilized 
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for observation of the data sets consisted of two software packages: The HYDICE Starter 
Kit, (an extension of) the Spectral Image Processing System (SIPS) and the Environment 
for Visualizing Images (ENVI). Both SIPS and ENVI provide graphical interfaces that 
utilize the Interactive Data Language (IDL) to perform spectral image processing 
functions. The HYDICE Starter Kit incorporates advanced detection algorithms not 
currently found in SIPS or ENVI. Similar equipment and software were used in the 
analyses done at The Naval Postgraduate School. 

a. HYDICE Starter Kit 

The author found that the HYDICE Starter Kit was extremely sensitive to the 
type of inputs or commands and many times would abort to the IDL prompt. This 
required exiting from IDL and rebooting the software each time. This could have been 
attributed to the author’s unfamiliarity with software operation, however, other members 
of the exploitation team have encountered the same problems. The SIPS View module 
of the Starter Kit can be utilized to visualize and perform limited analysis on large 
imaging spectrometer data sets. Color composite images are relatively easy to produce 
and extraction of individual spectra from pixels or selected groups of pixels is 
straightforward. Extracted spectra can be compared to previously collected laboratory or 
field spectra stored in a spectral library. The spectral library function allows the user to 
display and manipulate single/multiple spectra and to save the results for output. 

The Starter Kit incorporates advanced detection algorithms not currently carried 
as a part of SIPS, these are: Orthogonal Subspace Projection, Low Probability of 
Detection and Constrained Energy Minimization along with Principle Components. 

b. ENVI 

The author found the ENVI software package very easy to use and more 
forgiving of anomalous inputs. Image file information is easily displayed and ENVI 
allows multiple windows to be opened, thereby facilitating the comparison of selected 
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images. Plotting and annotation of spectra is straightforward and the software provides 
the capability to save the plots in a variety of output formats. 

2. Hyperspectral Image Data 

Hyperspectral data was collected over two areas. The first was the experiment area 
within the White Sands Missile Range complex and the second was over the urban and 
agricultural areas of Las Cruces, New Mexico. The objective of the data analysis was to 
determine if the detection algorithms would provide reliable indications of a target's 
presence, in this case, the VI tank underlying the F9 radar scattering camouflage. The 
collections over the Las Cruces urban area were not examined because locations of 
potential military related targets could not be confirmed and subsequently compared to 
detection images. In order to perform these analyses the experiment environment should 
be as controlled as possible. With this in mind, initial observations of the data sets were 
done to examine various aspects of the images. 

♦ The spatial quality of the images was sufficient to resolve features 

within the field of view. 

♦ The target in question was included in the image. 

♦ The inclusion of calibration targets within the same field of view as the 

target in question. 

♦ The spectral quality of the image. 

♦ The target/sensor geometry. 

The optical target described in Chapter IV was placed in the experiment area of Site 
B on 27 October, 1994. Each image was examined to see if the area in which the target 
was placed was in fact imaged by the sensors. 
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i. AAHIS Data 


Each hyperspectral data set was loaded and examined utilizing the HYDICE 
Starter Kit, for the above characteristics. In total, 29 sets of data were examined. Only 
one data set was found that met the above criteria. This set, the 27 October, 1994, Data 
set #8, was collected over site B at 2000" (AGL) and approximately 1250:51 hrs (MST). 
As shown in Figure 4.1, this scene includes the VI Tank target as well as the calibration 
panels. 

One data set was found that could possibly contain the optical target. The 
target's exact location was not known, however, an approximate location was determined 
after discussion with the experiment coordinator. The 28 October 1994 data set #13, as 
shown in Figure 4.2, showed potential for covering the target location. Spectral 
examination of the scene however showed that the AAHIS detector became saturated 
between approximately band 33 (.6129 microns) and band 53 (.7236 microns). This 
saturation made the image unusable for this analysis. 

b. GER (DAIS) Data 

The GER data sets were examined using the HYDICE Starter Kit. The images 
over the Site B area, collected on 28 October 1994 (data sets 11-14), were found to be 
severely undersampled. The undersampling, due to low sensor frame rates in relation to 
sensor platform velocity over the ground. As seen in Figure 4.3, this causes the image to 
take on an appearance that the scene is squashed and severely restricts any spatial 
resolution and hence, identification of potential targets. However, this fact does not 
effect the spectral quality of the image itself. It was possible to define some relatively 
linear features such as roads provided they were vertically presented in the image. These 
images were not considered in the analysis of the VI tank spectra because the spatial 
location of the tank target could not be identified. 

The optical target may have been sampled in the GER images, the problem was 
that the extremely poor spatial coverage of the site would not allow for verification of 
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potential detections in the general area of the target location. In addition, any false 
alarms, or false detections, would be difficult to isolate. 
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V. DATA ANALYSIS 


The objective of the analysis was to determine if an exposed target could be 
detected by using the Low Probability of Detection (LPD) or Constrained Energy 
Minimization (CEM) algorithms. In addition, a Principle Components (PC) 
transformation was conducted on the entire scene to examine any similarities that may 
exist between the LPD component image and any of the PC images. As described in 
Chapter IV, the 27 October 1994 AAHIS data set 8 was the only scene with the required 
spatial and spectral qualities to be utilized in the analysis. The specific object of interest 
was the M-60 tank, cataloged as vehicle VI, in Table 3.3. At the Site B set-up for 27 
October the VI tank target was covered with a radar scattering, desert camouflage net 
cataloged as F9, in Table 3.1. 

The analysis focused on whether the VI tank spectra could be detected through the 
camouflage. Since the openings in the camouflage were much less than the pixel size, the 
VI spectra would most probably be subpixel in nature and be mixed with the camouflage 
spectra in addition to any contaminant spectra from the vehicle. Note that in all the 
figures presented here, wavelength values are in microns and reflectance values are in 
percent reflectance. 

A. SPATIAL ANALYSIS OF THE SCENE 

The 27 October, data 8 scene is shown in Figure 4.1. The size of the scene is 192 
samples x 512 lines x 70 wavelength bands. For the scene presented, information for 
band 6 was arbitrarily selected (0.46207 microns). 

The general layout of the target area can be seen. Spatial resolution computed for 
the scene based on the 2000' (609 m) flight altitude and an IFOV of approximately 1.0 
milliradians, is 0.67 meters 2 . Visual quality of the image is good although there is some 
blurring present. The target area is bordered by the dirt road that cuts through the top 
right half of the scene and what looks like a burm located in the lower left portion of the 
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scene. The large (20' x 20') calibration panels are visible. The remaining target panels 
are located to the left of the calibration panels. Exact locations can be verified with the 
Site B target layout presented in Figure 3.3. 

The VI tank which is covered by the F9 camouflage is clearly visible in a blow-up 
of the scene as shown in Figure 5.1a. Figure 5.1b will be discussed later in this chapter. 
Heavy shadows are present on the north side of the tank indicating a midday sun position 
at collection time. However, it is not clear whether the shadow is caused by the tank 
body which is approximately 3 meters in height by 7 meters long or the covering 
camouflage. The tank presents a surface area of approximately 24 meters 2 on the top 
surface and 20 meters 2 along the north side while the camouflage has a total surface area 
of approximately 104 meters 2 . Very prominent tank tracks are visible along a path from 
the tank's location to the dirt road. Vegetation is sparse although there are a large number 
of Mesquite bushes located throughout the scene. 

B. SPECTRAL ANALYSIS 

The 27 October 1994, data 8 scene was loaded into the Hydice Starter Kit for initial 
analysis. The approach taken for the analysis was to assume that the nature and make-up 
of the prevailing background signatures present within the scene were unknown. In 
addition, the objective of the analysis was to detect a man-made target signature (VI tank 
spectra) that was situated within a naturally occurring background. Since the area of the 
exposed surfaces was small in relation to the overall scene, the VI tank spectra would 
most probably occur in relatively few pixels within the scene and not be statistically 
significant. For these reasons, the LPD algorithm was chosen for it's robustness in 
detecting sub-pixel targets given the above conditions. The CEM algorithm was not 
applied to the data because the VI Tank spectra did not constitute a sub-pixel target that 
occurred in a relatively large number of pixels within the scene. 



1. Data Conversion 


The ground truth data provided by TEC was sampled from 400-2500 nanometers 
with values represented in percent reflectance. Each target spectra was then resampled 
utilizing the AAHIS wavelength file to produce spectral plots of the target spectra within 
the AAHIS sensor range. 

Two files were created in ENVI utilizing the Empirical Line Method described by 
Farrand et al, (1994). The first was the gains file which incorporates the atmospheric 
transmittance and instrument factors and the second was the offsets file which takes into 
account the atmospheric path radiance. Both files were applied to the target spectra in 
order to convert them from surface leaving reflectance percentages to AAHIS sensor 
perceived radiance. Figures 5.2 and 5.3 show the gains and offsets plots. Both files were 
derived from samples taken from the calibration panels located at the middle right hand 
portion of Figure 5.1. Both files were constructed by Dr. William Farrand, a member of 
the HYDICE exploitation team. 

2. Data Preparation 

The methods described in paragraph B1 were applied to the selected target spectra. 
The next objective was to prepare the spectra for use with the LPD algorithm. The 
algorithm assumes that sensor noise and the combined at-sensor radiance of unmodeled 
components in the scene are Gaussian in nature, and are independent and identically 
distributed (iid). (Farrand and Harsanyi, 1995). When a data cube is converted to surface 
reflectance by application of gains and offsets files, the noise statistics of the image will 
be effected and the assumption that the noise is iid would not be true. For this reason, the 
target spectral signatures and the pixel values of the image should be represented in 
sensor perceived radiance. 

The ground truth spectra were converted to sensor perceived radiance by utilizing 
the spectral math capability of the Starter Kit. To convert each spectra, each ground truth 
data file was multiplied by the gains file and the product added to the offsets file. Figures 
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5.4 through 5.9 depict the original target spectral plots and the new target spectral plots 
after conversion to sensor perceived radiance utilizing the gains and offsets files. One 
should note the distinct similarities in the overall shapes of the spectral plots and the fact 
that the TANK spectra is an exact match to the VI TAN spectra. Figure 5.10 depicts the 
Site B soil spectra which can be used to compare to the target spectra similarities. 
Similarities in overall spectral shape can be seen out to 0.6 microns but the soil spectra 
takes on a more linear slope beyond that point out to the end of the sensor range. The 
target spectra all show distinct differences in the slope of their plots beyond the 0.6 
micron wavelength range and therefore should be discernible from the background. 
Figure 5.11 depicts the F9 Camouflage spectra for later comparison. 

3. Target Spectra Selection 

The LPD algorithm requires a target spectrum that is run against the entire data 
cube. TEC sampled five areas on the VI tank. The samples taken represented the brown, 
light tan, tan-brown painted areas of the tank along with the gray painted gun boot and a 
rusted area on the left rear of the tank. The plots of these various paint spectra along with 
an average of all the spectra collected over the tank is shown in Figure 5.12. Out of these 
spectral plots, three were chosen as representative of the VI tank spectra, these were: 

♦ The average VI spectra sampled and represented as TANK 

♦ The average VI Tan spectra sampled and represented as VI TAN 

♦ The average VI Tan Brown spectra sampled and represented as VI 
TAN BROWN 

To better compare the target spectra, Figure 5.13 has been provided. Due to the 
obvious similarity between the TANK spectra and the V1 TAN spectra, it was decided to 
use only the TANK and the VI TAN BROWN spectra as target signatures when running 
LPD. 
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Figure 5.14 compares the spectral shapes of the TANK spectra and the F9 spectra. 
Spectral similarities appear out to approximately 0.74 microns where the two plots begin 
to diverge from each other. These characteristic differences should allow for 
discrimination of the TANK spectra from the F9 camouflage. 

4. LPD Analysis 

The LPD algorithm was run with the TANK and the VI TAN BROWN spectra as 
target signatures. After selection of the scene and the target spectra, primary eigenvectors 
are computed and a correlation matrix constructed. In the context of the LPD algorithm, 
the primary eigenvectors characterize the majority of the spectral variability in the image 
and are used to form an operator that suppresses the undesirable background signatures. 
In the initial LPD runs, the first 11 primary eigenvectors were removed by the algorithm. 
It should be noted that if a certain material is present in an amount greater than the target 
spectra, but in a small enough quantity to be excluded among the primary eigenvectors, 
then it may show up as a false alarm in the LPD component image. One solution is to 
increase the number of eigenvectors to be removed, thereby increasing the amount of 
background that is eliminated. The LPD component image provides an assessment of the 
relative abundance of a particular material under the assumption that target and 
background spectra are linearly mixed. 

a. LPD Results (TANK Target) 

Figure 5.16 shows the results of the LPD run using the TANK spectra as the 
target signature. For this image, a 50% contrast stretch was applied. Figure 5.17 is the 
same LPD component image but with a histogram equalization applied. This is a special 
case of a contrast stretch, which enhances the detectability of the features within the 
image by mapping data values into display values such that a roughly uniform 
distribution of gray levels is presented. The VI TANK spectra shows up as a bright spot 
in the upper left comer of the image. This position can be compared to the location of VI 
Tank in Figure 5.1. Several false alarms occur within the image and are circled for 
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clarity. The brightest false alarm pixels indicating higher DN value assignment are 
evident in Figure 5.16 and come from the panels marked P3 (TAN CARC PAINT) and 
P4 (TAN TAGGANT). False hits also occurred on other panels such as the F2 (GREEN 
NYLON TENTING) and the F6 (TAN TWILL) and are more evident in Figure 5.17 after 
application of a Histogram equalization. This is most probably due to these spectra being 
included in the residuals but not included among the primary eigenvectors. Figure 5.18 is 
the same target signature run using 20 eigenvectors in order to eliminate more of the 
undesired background. The false alarms are more evident in this image due to contrast 
stretching that was applied to highlight the higher DN values. Overall, the LPD operator 
has performed well in suppressing the majority of the background in each of the 
component images. 

b. LPD Results (TANBROWN Target) 

Figure 5.19 shows the results of the LPD run utilizing the TAN BROWN target 
spectra and 11 eigenvectors. The differences between this component image and the one 
shown in Figure 5.16 are not readily distinguishable. This fact holds true again when 
comparing Figure 5.20, the same LPD image except with a histogram equalization 
applied, with Figure 5.17. Figure 5.21 is the LPD component image run for this target 
spectra using 25 eigenvectors. In this image, the abundance of the target spectra shown 
in the location of F9/V1 appears less. In addition, some of the false alarms on certain 
panels disappear while others appear for the first time. This may be caused by the 
elimination of certain residuals by removing more eigenvectors and the utilization of a 
different target spectra. A histogram stretch is applied again to this image for uniformity. 
Again, the LPD operator did a good job of background suppression using 25 
eigenvectors. 

Figure 5.15 compares four different spectra that show up as false alarms in the 
LPD runs. The P3 (TAN CARC PAINT) spectra shows the most spectral variability 
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across the sensor range, while the remaining three spectra, P9, F6 and Desert Color BDU 
(TAN) exhibit only minor slope/shape differences. 

5. Principle Components Analysis 

The Principle Components Transformation produces a set of uncorrelated images 
equal to the number of spectral bands in the image. The images with the most 
information or largest variance make up the first Principle Component with each 
successive image containing less variance. In addition, each successive Principle 
Component image will be orthogonal to the previous one. Typically, the first 6 to 10 
Principle Components will account for more than 99% of the total variance within a 
scene. Hyperspectral data, by nature, is highly correlated between each spectral band. 
This is due primarily to the fact that the spectrum of most materials vary slowly and in 
similar fashion. The main benefit in the use of the Principle Components Transformation 
with hyperspectral data sets has been in the area of data volume reduction. The approach 
taken for this analysis was to see if any of the Principle Component Images would closely 
resemble any LPD results and if so, what do the different components represent? 

A Principle Component Transformation was run on the 27 October data 8 set to see 
if any similarities exist between any of the principle Component band images and the 
LPD Component images. 

a. Principle Component Results 

The Principle Components Transformation was conducted using the ENVI 
software package. Each component image was examined for similarities with the LPD 
results. Figures 5.22 and 5.23 show component images 1 and 2. Component 1, 
containing the most variance within the scene, looks very similar to Figure 4.1. Figure 
5. lb is a blow-up of the VI Tank target from Principle Component 1. It is clear that the 
Tank target representation in this figure is more detailed than that of Figure 5.1a. 
Principle component 2 is distinctly different but bears no resemblance to either of the 
LPD component images. Component 8 and 10 showed the greatest likeness to the LPD 
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component images. Component 9 did not show any meaningful data. Figure 5.24 shows 
component 8. The F9/V1 tank target and panels F2, P3, P4 and P9 show up clearly in 
the image. No direct correlation with either LPD component image was observed. The 
10th component image shown in Figure 5.25, shows numerous target panels in addition 
to the F9/V1 tank target. A comparison with Figure 5.21, the LPD run using 25 
eigenvalues, reveals a close match on all target hits. This component should represent 
less than 1% of the total scene variance. This would indicate that the target panels are 
spectrally similar in this range and comprise a very small portion of the total pixels in the 
image. 
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VI. CONCLUSIONS 


Operation Desert Radiance was the preliminary attempt by the HYDICE Program 
Office in demonstrating the utility of hyperspectral imagers for detecting a tactical target 
based on it's spectral characteristics alone. The author's participation in the collection 
experiment and analysis of the sensor data has highlighted the complexities associated 
with an experiment of this size. 

The experiment was well designed but very ambitious. The use of the AAHIS 
sensor proved to be a good choice. However, the GER (DAIS) should not have been 
utilized based on it's poor signal-to-noise ratio and low frame rates. The use of "budget" 
pilots by the sensor operators severely reduced the usefulness of the data collected. 

Ground truth measurements during the experiment were performed by TEC and 
EG&G. The number of ground truth samples taken for each material appeared to be 
small. In addition, the samples were not well distributed over the entire surface of the 
tank and represented spectra of portions of the surface with various paint schemes. These 
spectra were averaged and may not be an accurate representation of the overall spectra 
from the tank. A more accurate method of collecting the ground truth would be to collect 
the entire spectra of the tank from an overhead position. It is understood, however, that 
the large number of materials to be sampled in such a short timeframe may have 
contributed to the small sample sizes. 

Distortions appeared in some of the calibration panels sampled by EG&G and the 
parachute target materials were not homogeneous in color as they were folded into panels 
on the ground. In addition, some warping of the painted panels was also observed. These 
facts, coupled with the fact that the majority of the samples were taken over the comers of 
the target panels, probably resulted the collection of spectra that were not representative 
of the entire panel. 

The results of the data analysis clearly show that unique target detection is possible 
utilizing the LPD algorithm. A higher than expected number of false alarms occurred 
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even though a large majority of the background was suppressed. The similarity between 
the Principle Component 10 image and the 25 eigenvector, LPD component image was 
intriguing. Both results show that the majority of the targets contribute very little to the 
overall scene variance and are subsequently seen only in the higher order component. 

The LPD algorithm performed well during the analysis but it's capabilities to detect 
a specific target spectra were limited because it could not utilize the spectral variability 
associated with wavelengths beyond 1.0 micron. The use of the HYDICE sensor in 
future collections should correct this problem by imaging out to 2.5 microns. 

The analysis shows that image processing algorithms for automatic target detection 
need more refinement to eliminate the high false alarm rate. Further target detection 
studies utilizing sensors that provide a greater degree of spectral resolution and which 
cover a broader spectral region will determine if the problems encountered in uniquely 
detecting a target are related to the spectral variability of the target material utilized in the 
algorithms. 
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APPENDIX. FIGURES 




WAVELENGTH (Nanometers) 


Figure 2.1 (Rinker, 1990) 
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Figure 2.5 (Kruse, 1987) 
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Figure 2.13 (Melack, Pilorz, 1990) 
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Figure 2.16 (Lehmann,et al, 1990) 
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Figure 2.22 (Cracknell and Hayes, 1991) 
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Figure 2.23 (Harsanyi and Chang, 1994) 



Figure 2.24 (Harsanyi and Chang, 1994) 



Figure 2.25 (Harsanyi and Chang, 1994) 
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Figure 4.2 (AAHIS 28 October Data 13) 
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Figure 5.4 (Tank spectra from TEC ground truth —AAHIS sensor range) 



Figure 5.5 (Tank spectra with AAHIS Gains /Offsets applied) 
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Figure 5.16 (LPD Component image-TANK Target Spectra) 
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Figure 5.17 (LPD Component image-TANK Target Spectra-Histogram Eq.) 
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Figure 5.18 (LPD Component Image-TANK I arget Spectra-20 Eigenvalues) 
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Figure 5.19 (LPD Component Image-TAN BROWN Target spectra) 










Figure 5.20 (LPD Component Image-TAN BROWN Target Spectra-Flistogram Eq.) 
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Figure 5.21 (LPD Component Image-TAN BROWN Target Spectra--25 Eigenvalues) 
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Figure 5 22 (Principle Component 1) 
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Figure 5.24 (Principle Component 8) 
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Figure 5.25 (Principle Component 10) 
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